Pis'ma V ZhETF 

Quantum electrodynamics with anisotropic scaling: Heisenberg-Euler 
action and Schwinger pair production in the bilayer graphene 

M. I. Katsnelson* i) and G.E. Volovik +* 

* Radboud University Nijinegen, Institute for Molecules and Materials, Heyndaalseweg 135, NL-6525AJ Nijmegen, The Netherlands 

+ Low Temperature Laboratory, Aalto University, School o{ Science and Technology, P.O. Box 15100, FL00076 AALTO, Finland 

* Landau Institute for Theoretical Physics RAS, Kosygina 2, 119334 Moscow, Russia 



(N 

o 

(N 



a^ 






I 

o 
o 



> 

00 

rn 

O 



X 



We discuss quantum electrodynamics emerging in the vacua with anisotropic scaling. Systems with 
anisotropic scaling were suggested by Hoi'ava in relation to the quantum theory of gravity. In such vacua 
the space and time are not equivalent, and moreover they obey different scaling laws, called the anisotropic 
scaling. Such anisotropic scaling takes place for fermions in bilayer graphene, where if one neglects the trig- 
onal warping effects the massless Dirac fermions have quadratic dispersion. This results in the anisotropic 
quantum electrodynamics, in which electric and magnetic fields obey different scaling laws. Here we discuss 
the Heisenberg-Euler action and Schwinger pair production in such anisotropic QED. 
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1. INTRODUCTION 

Both superfluid '^He-A [T] and single-layer graphene 
PJ 121 H] serve as examples of emerging relativistic quan- 
tum field theory in 3+1 and 2-f 1 dimensions, respec- 
'tively. Both systems contain "relativistic" fermions, 
iwhich are protected by the combined action of sym- 
metry and topology (see review 0), while the collective 
modes and deformations provide the effective gauge and 
.gravity fields acting on these fermions [I] [^ . 

Effective electromagnetic field emerging in super- 
fluid '^He-A is the collective field which comes from the 
^degrees of freedom corresponding to the shift of the po- 
sition of the Weyl (conical) point in momentum space. 
The corresponding Maxwell action (for special orien- 
tation of the effective electric and magnetic fields, see 
details in |T]) is 



S 



(fxdt 



Vp 

247r2 



B^- 



1 



-E' 



In 



1 



B2 



-^2 



(1) 

.where vp is the Fermi velocity. This action describes 
the effect of vacuum polarization, and corresponds to 
the Heisenberg-Euler action of standard quantum elec- 
trodynamics [B]. The only difference from the latter is 
that the fermions in the vacuum of "^He-A are mass- 
less, and their masslessness is protected by momentum- 
space topology of the Weyl point. This results in the 
logarithmic term which describes the zero-charge effect 
- the screening of the effective charge. The effective 
action becomes imaginary when E > vpB describing 
Schwinger pair production in electric field [7] with the 
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production rate oc E'^ for the case of massless fermions 
8^ (Schwinger pair production of massive fermions has 
been discussed for the other topological superfiuid - 
■^He-B, see Ref. [5]). Similar action, but for the real 
electromagnetic field, should take place in the Weyl 
semimetals (for semimetals with the topologically pro- 
tected Weyl points see Refs. [TUl HH [II [13 [H [IS] ) ■ 
Among the other things, the massless fermions give rise 
to the chiral anomaly, and the corresponding Adler-Bell- 
Jackiw equation for anomaly has been verified in the 
■^He-A experiments 16. . 

Graphene gives the opportunity not only to study 
the 2+1 quantum electrodynamics, but also extend the 
theory to different directions. In particular, to study 
systems with anisotropic scaling, which were suggested 
by Hofava for construction of the quantum theory of 
gravity, which does not suffer the ultraviolet (UV) di- 
vergencies (the UV completion of general relativity) 
[TTI ITSl [T9l [20] . As distinct from the relativistic mass- 
less fermions in a single layer graphene, which obey the 
invariance under conventional scaling r ^^ &r, t — > bt, 
fermions in a bilayer {N — 2) or rhombohedral A^-layer 
{N > 2) graphene have a smooth touching of the Dirac 
point, E^ ex p^^ , see Refs. [H[l]. These fermions obey 
the anisotropic scaling r — >■ &r, i — >■ b^t, precisely which 
is needed for construction of the divergence-free quan- 
tum gravity. 

Here we discuss the effect of this anisotropic scal- 
ing on the effective action for real or artificial (e.g., 
created by deformations [21 [1]) electromagnetic fields. 
Due to anisotropic scaling, which distinguishes between 
the space and time components, the electric and mag- 
netic fields experience different scaling laws. In par- 
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ticular, the one-loop action for the magnetic field is 
oc B^Pq ^ , where D is the space dimension and po 
is the infrared cut-off. This demonstrates that super- 
fluid ^He-A with its Z? = 3 and iV = 1 and the bi- 
layer graphene with its 13 = iV = 2 both correspond to 
the critical dimension Dc = 4: — N and thus they both 
give rise to the logarithmically divergent action for the 
magnetic field (which manifests itself as a logarithmic 
divergence of diamagnetic susceptibility in undoped bi- 
layer graphene with the trigonal warping effects being 
neglected [HI [23]). However, for the electric field the 
action is ex E'^p^^'^^^ . While for the superfluid ^He-A 
this again gives the logarithmic divergence as it happens 
in the relativistic systems with massless fermions, in the 
multiple-layered graphene with its D = 2 the action di- 
verges in the infrared for any A'^ > 1, giving rise to the 
power-law action ex i<;(^+2)/(-'^+i). 



2. EFFECTIVE ACTION FOR REAL AND 
INDUCED ELECTROMAGNETIC FIELD 

Anisotropic scaling for electromagnetic field 
The effective Hamiltonian for the bilayer graphene 
in the simplest approximation is [H H] 



in Ref. ^^. We consider the electrodynamics with 
anisotropic scaling, which is induced by fermions obey- 
ing the Hamiltonian ([5]). 

General effective action 

In Eq. ([2|) there is only one dimensional parameter, 
the mass m. Combining the dimensional analysis with 
the anisotropic scaling, one obtains that the effective 
action for the constant in space and time electromag- 
netic field, which is obtained by the integration over the 
2-1-1 fermions with quadratic dispersion, is the function 
of the dimensionless combination /x of electric and mag- 
netic field 
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where a are Pauli matrices and m is the mass en- 
tering the quadratic band touching. Experimentally, 
m Rs 0.03me where m^ is the free-electron mass ^^. 
Here we neglect the degrees of freedom related to the 
tetrad gravity and concentrate on the degrees of freedom 
corresponding to the electromagnetic field (in principle 
the vector potential A may include not only the real 
electromagnetic field, but also the collective field which 
come from the degrees of freedom corresponding to the 
shift of the position of the Dirac point in momentum 
space, as in "^He-A). We shall use the natural units in 
which ft = 1; electric charge e — I; the vector field 
A has dimension of momentum, [A] = [L]~^; electric 
and magnetic fields have dimensions [E] — [LT]^^ and 
[B] = [L]~'^, respectively. 

For standard quantum electrodynamics emerging in 
the vacuum with massive electrons, the Lorentz in- 
variance combined with the dimensional analysis gives 
the general form {B^ — E'^)f{x,y) for the Heisenberg- 
Euler action [5] in terms of dimensionless quantities 
x = {B^ - E^)/M^, y = B- E/M^. Here M is the rest 
energy of electron, which violates the scale invariance 
in the infrared. Extension of the Born-Infeld electrody- 
namics to the anisotropic scaling has been considered 
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The change of the regime from electric-like to magnetic- 
like behavior of the action occurs at /x « 1. The asymp- 
totical behavior in two limit cases, g{^ -^ 0) — a and 
g{li — > oo) — {b + ic)y?'^ , gives the effective actions for 
the constant in space and time magnetic and electric 
fields: 



SB = a I d^xdt — , Se^ I d^xdt{h + ic)E^/^m^/^ . 

(4) 

The dimensionless parameters a and h describe the vac- 
uum polarization, and the dimensionless parameter c 
describes the instability of the vacuum with respect to 
the Schwinger pair production in the electric field, which 
leads to the imaginary part of the action. 

This should be contrasted with the single layer 
graphene, where the corresponding effective action for 
the 2-|-l relativistic quantum electrodynamics is 
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critical field is Ec{B) = vpB; and the Schwinger pair 
production is ex E^/^ at B = (see Refs. [H [li [17] ) . 
Returning to the conventional units, one should substi- 
tute B -^ eB/c and E -^ eE. 

Action for magnetic field 

First, let us consider the case E = 0. The action is 



S 



d^xdt 



B^ 



(6) 



where a is undefined yet numerical coefficient. The sim- 
plest way to find a is the following. First, let us consider 
the case of uniform magnetic field and the sample of unit 
area, than, J d'^x —J' 1. Second, let us go to the imag- 
inary time, exp (iS) -^ exp (—j3F) where F is the free 
energy, /3 = 1/T is the inverse temperature. The change 
of the free energy in magnetic field is F = — x-B^/2 
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where x is the magnetic susceptibihty (per unit area). 

Thus, we have 

mX ,ys 

a=— (7) 

The susceptibihty for the case of bilayer within the 
model of purely parabolic spectrum has been calculated 
in Refs. [351 US]- As we know, the D = 2 vacuum of 
fermions with quadratic touching, N = 2, corresponds 
to the critical dimension for the magnetic field action: it 
is logarithmically divergent for the case of zero doping. 
This logarithm should be cut at the energy of the re- 
construction of the spectrum due to the farther hopping 
effects ("trigonal warping") and/or interelectron inter- 
action (for the most recent discussion, see Ref. [21] )■ 

The answer is 

9sgv 
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where gs = g^, = 2 are spin and valley degeneracies. 
For small fields, the logarithmically running coupling is 
A = In (7^/A^) , where the ultraviolet cut-off is provided 
by the interlayer hopping 71 , while the infrared cut-off is 
provided by the trigonal warping whose typical energy is 
A f« O.OI71. For larger fields, B > A^, the infrared cut- 
off is provided by the field itself, A = In (7^/5), which 
is similar to what occurs in effective electrodynamics in 
^He-A with massless fermions, see Eq.(IT]). 

Schwinger pair production 

Schwinger pair production in bilayer graphene in 
zero magnetic field has been considered in Ref. [27] . 
Let us consider the case of the crossed fields B = Bz, 
E = Eic using the semiclassical approximation. We will 
use the gauge A^ = 0,Ay = Bx. Thus, the imaginary 
part of the momentum along x direction is determined 
by the equation (cf. Ref. [H]): 



(x) = {k + Bx) - 2mE \x\ 



(9) 



where k is the (conserving) momentum in y-direction. 
The classically forbidden regions relevant for the tun- 
neling is determined by the condition k^ (x) > and 
the tunneling exponent is determined by the imaginary 
part of the action 



Se (fc) = 2 / dxK (x) 
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where Xl,Xj^ are the left and right turning points. 
One has to introduce the parameter 
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thus, 



where 
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Now we have to study conditions of existence of the 
left and right turning points. We can restrict ourselves 
by the case fc > only, due to symmetry with respect 
to the replacement k — >■ —k,x — >■ —x. Simple analysis 
show that two turning point exist only if 



£ >2k 



(15) 



The further calculations are straightforward, and the 
answer is 



Se (fc) == fJ-f 



J2k 
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The Taylor expansion for the function / reads 
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For f — i> 00 we have / (x) « x^/6, which is in agree- 
ment with Ref. [27], if one fixes the misprint in Eq.(36) 
of [22, where there is a power 1/3 instead of 3. The 
semiclassical approximation is valid if /i 3> 1 . 

The pair production is obtained after integration of 
the tunneling exponent over the momentum k: 
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Semiclassically, there is always some states which 
demonstrate tunneling, but only with |A:| < £/2. In the 
case of strong magnetic field their contribution shrinks 
to the point. 

In the limit of small magnetic fields one obtains the 
parameter c in the action ^: 



IniS — ' / ^ B^u. / dxexp { ——x^] 
Att^ Jo V 6 / 



ci?4/3mi/3 . 



c = M!_6i/3r^/3) (20) 
127r^ 
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The Schwinger pair production rate proportional to 
£;''/3 has been discussed in Ref. [27] . 

If we take into account the next-order corrections 
in 1//X the expression (pn|) is muhiphed by the factor 

1 25/322/3 j^2 

-*- ~ 10r(l/3) {mE)*/3- 

The change of the regime from magnetic-hke to 
electric-hke happens at /z « 1 which means, in CGSE 
units, 



n = (T+ {{ei + ie2) ■ {p - eA)Y 
+(T" ((ei - ze2) ■ (p - eA))^ 



(25) 
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For the field i? = IT this ratio is of the order of 3 • 10~^, 
that is, the order of magnitude smaller than for the case 
of the single-layer graphene where it is vp/c « 1/300. 
Effective aetion for higher order touching 
All this can be extended to the order of N band 
touching, which presumably may be achieved by rhom- 
bohedral stacking of N graphene- like layers [2j |4] . The 
effective Hamiltonian for fermions in this case is 

^ = ^ ((i + *y) • (P - ^A))^ + |- ((i - .y) . (p - eA)) 

(22) 
If we are interested only in the infrared behavior of the 
action in the regime when the electric field dominates, 
then the induced electromagnetic action has the follow- 
ing general structure: 



where ei and e2 play the role of zweibein fields, which 
give rise to the energy spectrum E"^ = {g^^PiPk) corre- 
sponding to the effective 2D metric g^^ — e\e\ + e\e\; 
to spin connection and torsion. 

As distinct from superfluid ^He-A, which is the 
analog of the relativistic vacuum with massless Weyl 
fermions, the bilayer graphene is the representative of 
the quantum vacua, which experience different scaling 
laws for space and time. While such vacua were con- 
sidered by Hof ava in relation to quantum gravity, here 
we applied the anisotropic scaling to quantum electro- 
dynamics emerging in these systems, using as an exam- 
ple the 2D system with massless Dirac fermions with 
quadratic spectrum. Such systems have peculiar prop- 
erties, and we touched here the Heisenberg-Euler action 
and the Schwinger pair production. 
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At large ^ one has g{fi) ex ^(^+2)/2(JV+i)^ and the 
Schwinger pair production rate in zero magnetic field 
is 

1 W + 2 

n - — (m^)«+i . (24) 

m 



3. DISCUSSION 

In the condensed matter context (or in related mi- 
croscopic theories of the quantum vacuum), the exis- 
tence of the topologically protected nodes in spectrum of 
Weyl media gives rise to the effective gauge fields (C/(l) 
and SU{2)) and gravity. The same takes place in single 
layer and bilayer graphene, which contain Dirac points 
in their spectrum. The effective SU{2) gauge field comes 
from the spin degrees of freedom, these are the collective 
modes in which the momentum of the Dirac point shifts 
differently for spin-up and spin-down species. In addi- 
tion to the collective modes related to the shift of the 
Dirac point, the graphene has degrees of freedom, which 
correspond to tetrads in anisotropic gravity. In partic- 
ular, for the bilayer graphene these degrees of freedom 
enter the effective Hamiltonian in the following way 
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